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ABSTRACT 

The transient, ultra-luminous X-ray source CXOM31 J004253. 1+41 1422 in the Andromeda 
galaxy is most likely a 10 solar mass black hole, with super-Eddington luminosity at its peak. 
The XMM-Newton spectra taken during the decline then track luminosities of 0.86 -0.27LEdd- 
These spectra are all dominated by a hot disc component, which roughly follows a constant 
inner radius track in luminosity and temperature as the source declines. However, at the high- 
est luminosity the disc structure should change due to advection of radiation through the disc. 
This advected flux can be partly released in the plunging region, with the remainder being ad- 
vected below the event horizon. We use a fully relativistic disc model, SLIMBH, that includes 
these effects, and uses full radiative transfer through the photosphere based on TLUSTY. It 
also incorporates full relativistic photon ray-tracing from the proper location of the disc pho- 
tosphere rather than the mid-plane as the disc is no longer geometrically thin. However, we 
find that these new models differ only very slightly from the non-advective (standard) BHSPEC 
even at the highest luminosities considered here. While both discs can fit the highest luminos- 
ity data, neither is a good fit to the lower luminosities. We speculate that this could be due to 
a decreasing fraction of magnetic pressure support with increasing luminosity, such as seen 
in models where the magnetic pressure saturates to some fraction of the gas rather than total 
(gas plus radiation) pressure. 
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1 INTRODUCTION 

The most luminous objects in the Universe are powered by accre- 
tion of matter onto compact objects like neutron stars and black 
holes. The energy released depends on the mass and spin of the 
central black hole as well as the mass accretion rate through the 
disc, but there is a maximum luminosity for which gravity is able 
to exceed the outward pressure of radiation. Ultra-luminous X- 
ray sources (ULX) are defined as objects whose bolometric lu- 
minosity exceeds this Eddington limit for a ~ 1OM black hole 
(~ 1.3 x l(fi 9 erg s~ [ ). These objects could be powered by nor- 
mal (sub-Eddington) accretion onto intermediate mass black holes 
(Mbh — 10 2 - 10 5 M o ). However, there are theoretical probl ems 
concerning the formation of such objects dKing et all 1200 lh . so 
these are unlikely to form the bulk of the ULX population. In- 
stead, the majority of sources are probably .super-Eddington ac- 
creting stellar mass black holes. Such flows come in two flavours, 
either exceeding the Eddington limit by power i ng strong out- 
flows { IShakura & Sunvaev|[l973l ; lLipunovd]l999l ; IPoutanen etafl 
1 2007th or by advecting the radiation along with the flow (Po l- 
ish doughnuts: I Abramowicz et al1dl978 ): Jarosz vnski et aljjl980h : 
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slim discs: lAbramowicz et al.l dl988b; ISad owski l2009h ). These 
processes can occur together l lPoutanen et alj2007T) . a s shown in re- 
cent numerical simulations of supe r-Eddington flows dOhsuga et al.l 
120091 : lOhsuga & Minesh"ige1l201 l|). In all these cases, radiation is 
emitted at the local Eddington limit, so that the total source lumi- 
nosity is L ~ Le<m(1 + lnAf/Mndd)- Advective accretion discs that 
locally can exceed the Eddington lim it and thus excite winds were 
proposed bv lDotan & Shaviv f20ll . 



The BHSPEC models are the best current su b-Eddington accre - 
tion flow models for stellar mass black holes l lDavis et"al]|2005t) . 
These calculate the spectrum from full radiative transfer through a 
disc atmosphere for a Novikov-Thorne (general relativistic) emis- 
sivity, including ray tracing incorporating general and special rela- 
tivistic effects. These were modified to include the effects of advec- 
tion bv lBursa et ail | |2012| . in preparation). Advection sweeps some 
of the emitted energy along with the flow as photons are trapped 
in the optically thick disc. However, these photons can be released 
again at lower radii as the material accelerates towards the black 
hole after the inner disc edge. Pressure forces can change the posi- 
tion of this inner edge from that expected from the innermost stable 
circular orbit, and increase the height of the photosphere as the disc 
is no longer flat. The resulting SLIMBH models were fit to spectra 
from the accreting black hole in LMC X-3. However, the differ- 
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ences between BHSPEC a nd SLIMBH are not l arge here although this 
object reaches ~ 0.6L Edd dstraubetalj|201lh . 

Here instead we fit these new models to the higher Eddington 
fraction flows seen in ULX. In particular, we use the XMM-Newton 
spectrum from the transient ULX CXOM31 J004253. 1+41 1422 
(hereafter M31 ULX-1). The advantage of this object is that it 
showed a steady exponential decline in luminosity in five XMM- 
Newton spectra from 1.6 - 0.8 x 10 39 erg s~ l after its discovery at 
5 x 10 39 erg s~ l in Chandra imaging data. Such exponential decays 
with this time scale are well known from transient low mass X-ray 
binary black holes in our Galaxy , making it most probable that this 
ULX is also a 10M Q black hole dMiddleton et alj|2012l ; lKaur et al.l 
l2012h . This mass would then imply that the luminosity sampled 
by the XMM-Newton spectra corresponds to 0.86 - O.HLem, con- 
siderably higher than seen in LMC X-3. Thus it provides an ideal 
testing ground for the slim disc models, especially as the Galactic 
absorption column is fairly low, allowing the broad band disc shape 
to be seen down to low energies. 

The paper is structured as follows. The origin and analysis of 
our data is specified in section(2] the modelling of the ULX spectra 
is described in section[3]and the results are discussed in section[4] 



2 X-RAY DATA: REDUCTION AND ANALYSIS 

To ensure consistency, we followe d the procedure of data reduction 
detailed in lMiddleton et al] d2012h and extracted the event files us- 
ing sas vlO and filtered the products for standard patterns (<=12 for 
MOS, <= 4 for PN) and flags (=0 for spectral and timing products). 
The full field hard (10-15 keV) count rate was used to create good 
time inter vals excluding contam ination by soft proton flares (see 
Table 1 o f lMiddletonetai]|2012h . These were then used to extract 
spectral (including responses) and timing products from a circular, 
35" radius source, with background regions taken from the same 
chip, avoiding other sources in the field. 

The first observation was slightly piled up (seen using the tool, 
epatplot), so the inner 5" centroid centred on the source was re- 
moved in the PN. The MOS is more affected by pileup, so was 
discarded for this observation. 



3 MODELLING THERMAL ULX SPECTRA 

For the trans-Eddington luminosity regime of ULX spectra the 
best available mo del to date is the relativistic slim disc, SLIMBH 
dBursa et al .120121 . in preparation). It accounts for effects connected 
to high mass accretion rates such as advection of radiation, the relo- 
cation of the innermost stable circular orbit (ISCO) towards smaller 
radii and the proper, non-equatorial location of the disc photosphere 
(see ISadowski et alj[2oTTI ; Istraub et alj|201 ll) . For the vertical ra- 
diative transfer SLIMBH di rectly incorporates the TLUSTY grid of 
disc local annuli spectra dHubenv & Land 1 19951) which are then 
integrated and ray-traced through the relativisti c spacetime. It is 
thus c onstructed in the same manner as BHSPEC dDavis & Hubenvl 
120061) . We note that the difference between the two models is 
the unde rlying accretion disc: Whi le BHSPEC is based on a stan- 
dard thin lNovikov & Thome! < fl973l) disc, SLIMBH uses a slim disc 
dSadowskill2009l) . A comparison of these two models offers valu- 
able insi£hJ_£nJte_effec^of advection (see section [XT}. 

iMiddleton et al.l j2012h fit the five XMM-Newton spectra of 
M31 ULX-1 simultaneously with TBABS x BHSPEC using the 
neutral hydrogen column density N H = 6.7 x 10 20 arT 2 (after 
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Figure 1. X-ray spectra of M31 ULX-1, modelled with SLIMBH. Spectrum 
and data/model ratio for the highest (= 0.85Lgdd. black) and the lowest 
(= 0.27LEdd, cyan) luminosity data are shown. 

Dickey & LockmarJ|l990l) . distance D = 780 kpc dVilardell et all 
Tanaka et al.ll2010h and viscosity parameter a = 0.01. With 
black hole mass fixed to the provisional value of M = 10M o 
they derived best fit inclination of i = 30° and dimensionless 
spin parameter a, = 0.36^JJ | ( j ) . The fit is fairly good, with xl = 
3169.4/2801, though the spectra are significantly better fit with a 
phenomenological model of TBABS X (DISKPBB+COMPTT) (xl = 
2752.2/2777). 

We now replace the BHSPEC model with SLIMBH, again with 
a = 0.01, to see if this can better reproduce the spectral shape, and 
give a physical model which fits the data as well as the phenomeno- 
logical one. Advection should become increasingly important after 
L ~ 0. lLedd so the differences between BHSPEC and SLIMBH could 
be significant. We find that the more physical models do fit the 
data slightly better, but that the improvement is only small, with 
xl = 3132.2/2801 for the same resulting black hole spin (see Table 

LB- 

A closer i nspec tion of the spectral residuals shown in 
IMiddleton et al.l d2012l) explains this result. The residuals between 
BHSPEC and the data are largest at lowest luminosities, with the 
data being increasingly well fit by the non-advective disc models 
at higher luminosities. This is opposite to the behaviour expected 
if advection becomes increasingly important, limiting the impact 
of the spectral changes from SLIMBH. We explore the differences 
between BHSPEC, SLIMBH and the data in more detail below. 



3.1 Effect of advection: The difference between SLIMBH and 
BHSPEC 

Figure [2^ shows the best fit model spectra of SLIMBH (red) 
and BHSPEC (blue). Their ratio (black/cyan) for identical mass 
(10M G ), spin (a, = 0.36), inclination (30°) and luminosity (= 
0.86/0.27L Ed d) is presented in the top panel of Figure[2p. For the 
highest luminosity, the slim disc is very slightly softer below the 
peak as radiation is transported inwards at these radii rather than 
being emitted. Some fraction of this advected radiation can then 
be released below the ISCO, so there is additional radiation in the 
SLIMBH models at the highest energies i.e. above the peak. Thus 
the SLIMBH models are broader and not so peaked as the corre- 
sponding BHSPEC models, but the difference is less than a few per 
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Table 1. Slim disc model: best fitting spectral parameters 



TBABS * SLIMBH 


Observation 1 


2 3 


4 5 


Nh [cnr 2 ] 
Mbh [M q ] 

c, 

a 


6.7 x 10 20 
10 

36 +0 02 

n cqq+0.006 n /ns+0.006 
°' 593 -0.005 °- 428 -0.006 

0.01 


n "«fi+0.004 a 971 +0.004 
-0.004 ' -0.004 



in 

D [kpc] 



30 
780 



x/dof<x$) 



3132.2/2801 (1.12) 



cent for 0.86L Ed(i . Advection only makes a more substantial differ- 
ence to the emitted spectrum for super-Eddington luminosities, and 
there is even less difference between SLIMBH and BHSPEC at lower 
luminosities (0.27LEdd)- 

However, there is a significant change in how well the models 
describe the data. The green lines in Figuref2ji show the best fit phe- 
nomenological DISKPBB+POWERLAW model for the highest/lowest 
luminosity dataset. The highest luminosity dataset is well described 
by either disc model. By contrast, the lower luminosity dataset is 
clearly broader than the corresponding BHSPEC and SLIMBH mod- 
els. The discrepancies between physical and phenomenological 
models are illustrated in the bottom panel of Figure |2p where the 
model ratio SLIMBH vs DISKPBB+POWERLAW for the highest (black) 
and lowest (cyan) luminosity is given. 

Thus, including adv ection does not change the conclusion 
of iMiddleton et alj J2012T) that the data (especially at lower lumi- 
nosity) are better fit by phenomenological two component models 
rather than our best current physical disc models. This shows that 
there is something missing from our best current disc models. We 
stress that the level is only 5-10 per cent, but this is significant given 
the quality of data now available. The missing physics is not ad- 
vection, nor is it any of the other spectral distortions which might 
become more important as the luminosity increases, such as winds, 
or bulk motion turbulence in the disc, as the data become better fit 
by our disc models as the luminosity approaches Eddington. 



4 DISCUSSION AND CONCLUSIONS 

It has long been assumed that in the trans-Eddington luminos- 
ity regime two effects could become i ncreasingly important for 
accretion disc model s, (i) advection jAbramowicz et al" 1988 



Mineshige et alj 2000l) and/or (ii) outflows llShakura & Sunvaevl 
19731 : IPoutanen et alj|2007l ; lOhsuga et al]|2009l) . In this paper we 
pin down the effect of advection by fitting ULX spectra with slim 
discs. We find that advection typically changes the best fitting spec- 
tral shape by removing photons in the intermediate energy range 
where the spectral peak is located and by releasing more photons 
at high energies (see in Figure|2p, upper panel). This is the charac- 
teristic signature of slim discs where photons are trapped in the ac- 
cretion flow, accreted inward and released from inside the plunging 



region (see I Abramowicz et al.ll20l"ct for a discussion of the inner 
edge of slim discs). 

However, these changes make very little difference to the emit- 
ted spectrum for L < L^aa, and we show here that this cannot 
be the reason for the mismatch between data and models seen in 
moderate luminosity (0.1 - 0.5LEdd) sp ectra from both this ULX 
and from galactic black hole binaries dKolehmainen et ail 1201 it 
IStraub & Ghasemnezhadll2012l) . Real accretion discs at these mod- 
erate luminosities have a spectrum which is subtly broader than 
is currently expected from the best disc models. This additional 
broadening is most marked at lower mass accretion rates and gets 
weaker as the disc luminosity increases. This is exactly opposite 
to the luminosity behaviour expected from either winds, advection, 
bulk turbulence and/or inhomogeneities in the disc. 

Instead, it could be connected to the (low level) coronal emis- 
sion. We cannot constrain a weak high energy tail in our data, 
but this could potentially illuminate the disc and hence change 
the structure of the photosphere. How ever, we note that the mos t 
disc dominated GX339-4 spectrum in iKolehmainen et ail d201 ll) . 
where the tail carries less than a few per cent of the total bolomet- 
ric flux, shows this additional broadening. This is a very low level 
of illumination to cause a noticeable change in disc structure. Self- 
illumination of the inner disc by light-bending around the black 
hole may be a larger effe ct when the coronal emission is low (e.g., 
iMiniutti & Fabianll2004f) . However, this should involve a constant 
fraction of the emission, rather than give a decreasing fraction as 
LI Leu increases. 

Here we speculate that this may actually be a signature of 
magnetic pressure support in the disc. This will increase the scale 
height of the disc over and above that expected from a standard 
disc, so decreasing its density and hence its true absorption opac- 
ity. Electron s cattering then bec omes more important, broadening 
the spectrum jDavis et"ai]|2009h . The extent of this change in disc 
structure depends on the ratio of magnetic pressure to total (radi- 
ation plus gas plus magnetic) pressure. This could explain the ob- 
served trend in the data if the magnetic pressure saturates to some 
fraction of the gas pressure rather than total pressure, as it would 
decrease in importance as L/Lem increases. 

Other observations also strongly indicate that there is a change 
in behaviour in magnetic stress scaling between the gas pres- 
sure and radiation pressure dominated regimes. The rapid rise 
to outburst, which occurs (mostly) in the gas pressure domi- 
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Figure 2. Left: The SLIMBH (red) vs BHSPEC (blue) models for a high/low (0.86/0.27Le ( j ( j) luminosity disc. The black hole mass, spin and inclination are all 
tied at 10M Q , 0.36 and 30°, respectively. Advection makes little difference (less than a few percent) at even the highest luminosity considered here. Significant 
changes to the emitted spectrum due to advection are only seen at super-Eddington luminosities. The best fit phenomenological models to the data (green) show 
that the data are well described by the physical disc models at the highest luminosity, but are significantly broader at lower luminosity. Right: The corresponding 
model ratios SLIMBH vs BHSPEC (upper panel) and SLIMBH vs DISKPBB+POWERLAW (lower panel) for the same high/low (black/cyan) luminosity disc. 



nated regime, requires very efficient transport of angular momen- 
tum as parametrise d by a Shakura-Sunyaev a viscosity of ~ 0.1 
jPubus et al.lll999h . However, if this efficient transport is main- 
tained into the radiation pressure dominated regime then the disc 
becomes effectively optically thin, and its colour temperature cor- 
rection increases markedly. This is in sharp contrast to the data, 
where the observed constancy of the colour temperature correction 
requires a < 0.01 at high L/L E dd , where the disc is in the radia- 
tion p ressure dominated regime jDone&DavIsll2008l ;l Straub et al .1 
1201 ll) . Thus the data from the ULX reported here, and previous 
work on the disc spectra seen from black hole binaries both sup- 
port a change in magnetic stress scaling with total pres sure. Current 
nume rical simulations do not show this behaviour dHirose et al.l 
120091) , but these are at the limit of computational capabilities. Fu- 
ture research on how the self generated magnetic field from the 
MRI scales with pressure could resolve this issue. 
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